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Correlation Matrices of Doubly Selective Fading 
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Abstract — In this paper, the Cramer-Rao lower bound (CRLB) 
of the sample frequency correlation matrices (SFCM) is derived 
based on a rigorous model of the doubly selective fading channel 
for orthogonal frequency division multiplexing (OFDM) systems 
with pilot-symbol-aided modulation. By assuming a fixed pilot 
sequence and independent samples, SFCM is complex Wishart 
distributed. Then, the maximum likelihood estimator (MLE) 
and the exact expression of CRLB are obtained. From CRLB, 
the lower bounds of total mean squared error (TMSE) and 
average mean squared error (AvgMSE) independent of the pilot 
sequence are deduced, which reveal that the amount of samples is 
the dominant factor affecting AvgMSE while the signal-to-noise 
ratio and the maximum Doppler spread have negligible effect. 
Numerical simulations demonstrate the analytic results. 

Index Terms — CRLB, Frequency correlation matrix, Doubly 
selective fading channels, OFDM. 

I. Introduction 

Playing a key role in the channel estimation for orthogonal 
frequency division multiplexing (OFDM) systems, the fre- 
quency correlation matrix (FCM) is utilized by many statistics- 
based channel estimation algorithms, e.g., the linear mini- 
mum mean-squared error (LMMSE) estimator and its optimal 
low-rank approximations [1], the MMSE estimator exploring 
both time and frequency correlations [2], the two-dimensional 
Wiener filtering [3], and those algorithms based on parametric 
channel model [4] [5] [6]. In real applications, the sample 
FCM (SFCM) is used in stead of the true one, and usually 
obtained through the least squared (LS) channel estimation. 

For fixed or slowly moving radio channels whose Doppler 
spreads are relatively small, the channels reveal a feature 
of block-fading [7]. Hence, the Doppler spread affects the 
accumulation of SFCM negligibly. However, for fast moving 
radio channels, intra-symbol fading becomes so prominent 
that inter-carrier interference (ICI) takes effect by not only 
degrading the performance of OFDM systems [8], but also 
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affecting the mean and covariance of SFCM. In [9], the bounds 
of the ICI power is derived. 

The bias-property of SFCM in doubly selective fading 
channels has been investigated in [10]. As a counterpart, in 
this paper, we find out the Cramer-Rao lower bound (CRLB) 
for FCM to evaluate the performance of maximum likelihood 
estimator (MLE) and uncover the factors influencing the 
estimation accuracy. 

This paper is organized as follows. In Section |II] the OFDM 
system and channel model are introduced. Then, in Section HIH 
CRLB for FCM is derived and further discussed to uncover 
the essential factors. Numerical results appear in Section |IV] 
Finally, Section [V] concludes the paper. 

Notation: Lowercase and uppercase boldface letters denote 
column vectors and matrices, respectively. (•)*, (•)^, and 
\\-\\f denote conjugate, conjugate transposition, and Frobenius 
norm, respectively. ^ denotes the Kronecker product. E{-) 
represents expectation. [A]i,j and [a]i denotes the {i,j)-th 
element of A and the i-th element of a, respectively, diag(a) 
is a diagonal matrix by placing a on the diagonal. 

II. System Model 

Consider an OFDM system with a bandwidth of BW = 
l/T Hz (T is the sampling period). N denotes the total number 
of tones, and a cyclic prefix (CP) of length Lcp is inserted 
before each symbol to eliminate inter-block interference. Thus 
the whole symbol duration is Tg — {N + Lcp)T. 

The complex baseband model of a linear time-variant mo- 
bile channel with L paths can be described by [11] 



(1) 
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where t; E TZ is the normalized non-sample-spaced delay 
of the l-th path, and hi{t) is the corresponding complex 
amplitude. According to the wide-sense stationary uncorre- 
cted scattering (WSSUS) assumption, /i;(t)'s are modeled as 
uncorrected narrowband complex Gaussian processes. 

Furthermore, by assuming the uniform scattering environ- 
ment introduced by Clarke [12], /i;(i)'s have the identical 



normalized time correlation function (TCF) for all I's, thus 
the TCF of the /'s path is 

rtA^t)=af Jo (27: fdAt) (2) 

where af is the power of the l-th path, is the maximum 
Doppler spread, and Jo( ) is the zeroth order Bessel function 
of the first kind. Additionally we assume the power of channel 
is normalized, i.e., X^^^o^ '^f = ^■ 

Assuming a sufficient CP, i.e., Lcp > L, the discrete signal 
model in the frequency domain is written as 

y/ (") = H/ (n)x/ (n) + Uf (n) (3) 

where x/(n), y/(n), n/(n) G C^^^ are the n-th transmit- 
ted and received signal and additive white Gaussian noise 
(AWGN) vectors, respectively, and H/(m) e qNxn 
channel transfer matrix with the (fc + i^, fc)-th element as 

JV-l L 
m=0 /=1 

where hi{n, m) — hi{nTs + {Lcp + m)T) is the sampled com- 
plex amplitude of the l-th path, k and w denote frequency and 
Doppler spread, respectively. Apparently, as H/(n) is non- 
diagonal, ICI is present. In fact, when the normalized max- 
imum Doppler spread fdTs < 0.1, the signal-to-interference 
ratio (SIR) is over 17.8 dB [13]. 

III. CRLB FOR Frequency Correlation Matrices 

Usually SFCM is obtained through the LS channel esti- 
mation. We consider OFDM systems adopting pilot-symbol- 
assisted modulation (PS AM) [1], hence only pilot symbols, 
denoted as yp{n) G C^^^, are extracted and used to perform 
LS channel estimation. In addition, the pilot sequence is 
assumed to be invariant along the time. Therefore, 

h.p,is{n) = Xp Vp(") = Xj;^Hp(n)xp + ^np(n) (5) 

where Xp = diag(xp) is a diagonal matrix consisting of pilot 
symbols, and the noise term is np{n) ~ CA/'(0, cr^Ijv)- 

Denote the instantaneous channel impulse response (CIR) 
vector as ht(n,m) — [hi{n,ni), . . . , hL{n,ni)]'^ , m = 
0, . . . , N ~ 1, according to the assumptions of WSSUS and 
uniform scattering, h.t{n,m) is complex normal, i.e.. 



where D = diag(cr;^), / = 1, , 



ht{n,m) -C7Vl(0,D) 

Tf),l — 1,...,L. Then form the CIR matrix 
as Ht(ri) = [hi(n, 0), . . . , ht{n, N — 1)], so we have 

vec(Hf(n)) - CA/'ljv(0, O ® D) 

where O G qNxn ^ Toeplitz time correlation matrix (TCM), 
defined as 



\ft]. 



= Jo (27r/d(mi - m2)T) 



(6) 



Then according to (|4]i, the channel transfer matrix H/(n) = 
F7-Hf(n), where G C^^^ is the unbalanced Fourier 
transform matrix, defined as [Fi-lfc.; — g^-'^'^fe'^!/^. Thus 



H/(n) ^ CA/'ArxAr(0,f2 ® (F.DFf )) 



Assuming CIR is independent of the thermal noise, with Q 
and (|7]i, we have 



hp,i,(n) ~CA/'jv(0,S) 
where the covariance matrix S is defined as 



s = c.Xp-i(Rp + ^i^)x;^ 



(8) 



(9) 



where cu = x^f Oxp, and Rp = Ft-DF:^ is the true FCM. 

When the LS estimated CFR's, i.e., hp ;s(n)'s, are available, 
SFCM is constructed as 
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(10) 



where Nt is the amount of samples. To derive the probability 
density function (PDF) of SFCM, we assume that samples are 
independent of each other, which may be a strict constraint. 
However, when the maximum Doppler spread is large and the 
spacing between two contiguous pilot symbols is compara- 
tively small, the correlation between them is rather low, which 
alleviates the effect of model mismatch. Then, based on the 
assumption of independence and ([8]), we know that SFCM has 
the complex central Wishart distribution with Nt degrees of 
freedom and covariance matrix S' = H/Nt [14], denoted as 



R 



■p,ls 



'CWjv(iVt,S') 



and its PDF is 

/(Rp,is) ^ 



etr(-S^-iRp,;,)(det(Rp,;,))^'-^ 



(11) 



(12) 



Cr^(iVt)(det(S')^* 

where etr(-) — exp(tr(-)) and Crjv(iVi) is the complex 
multivariate gamma function, defined as 

N 

cr^{Nt) = II r{Nt -k + i) 

k=l 

Then, from ( fT2] i. the likelihood function is written as 

L(Rp) = tr(-S'-iRp,,,) + {Nt - TV) ln(det(Rp,,,)) 

- HCTNiNt)) - Nt ln(det(S')) 

Therefore, the score function with respect to the parameter 
matrix Rp is 

scorefR ) - - ^^^^ X (13) 

score^Hpj - g^^^^j^^^ - ^yg^lRp) 9vec(S') ^ ^ 

where the first term on the right-hand side of ( fT3] l is 

9vec(S0^ _ ^ 
avec(Rp) ~ TVt 
and the second term is 

. vec[(S'-Rp,.S'- - Nti:'-r] (15) 

By letting the score function equal zero and with (|9|, the 
MLE of FCM is derived as 



-(Xp-^®Xp-i) 



(14) 



(7) 



MLE(Rp) 



XpRp^fsX^ — (j'^In 



X^l^Xp 



(16) 
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Fig. 1. Fitting tlie maximum eigenvalues of f2 witli {21\ for different sizes 
of FFT (N) and normalized Doppler's (JdTs). 



Note that ( fT6b relies on the pre-known TCM, i.e., fJ, and noise 
power. For Rayleigh fading channels, it means the maximum 
Doppler spread, fd, is known. 

Further, according to the score function, the Fisher Infor- 
mation matrix with respect to Rp is 



J(Rp) = E 



( 9L(Rp) \ / aL(Rp) 



\9vec(Rp) / \9vec(Rp) 

With OCllKinil, we have 
9L(Rp) 



H 



(17) 



9vec(Rp) 



xvec[(S'-iRpj,S'-i 



SO, ( [TT] ) is rewritten into (fTSb . shown at the bottom of the next 
page. Notice that 



and 



therefore 



i?{vec[(S'-iRp,,S'-i - 7VtS'-i)^] 

xvec[(S'-iRpj,S'-i - NtT.'-^f]"} 
= Var{vec[(S'-iRp,i,S'-i)^]} (19) 

Given S ~ CWjv(^t, S'), the entry of its second origin 
moment is [15] 



Therefore, the entry of its second central moment is 

i;[([s],,, -ii;([s],,,))([s]M -^^([s]m))] = 

Accordingly, ( fT9l l is rewritten as 

Var{vec[(S'^iRp,i,S'~i)^]} = iVt(S'-^ ® S'-^) (20) 



Then, with (|201l, J(Rp) is 



j(Rp) = -(x^-^®x;i)(s'-^ 



)(Xpi®x;^) 



From the Fisher Information matrix, the CRLB of Rp can 
be derived as [16] [17] 

CRLB(Rp) = J^^Rp) 

= ^(Xp ® Xf ® S'^)(Xf ® Xp) 

= -^(iXpS«Xf)®(ixfS^Xp) (21) 
With (|9]l, (ISTT l can be further written as 



CRLB(Rp) 



1 2 2 

— (Rp + ^lAr)®(Rp + ^I^f (22) 



Based on (|22| |. a lower bound of the total mean squared 
error (TMSE) for MLE(Rp) is 



TMSElb(Rp) = ti-(CRLB(Rp)) 



(23) 



where 7 = cr~^ is the signal-to-noise ratio (SNR). And, 
accordingly, the lower bound of the average mean squared 
error (avgMSE) is 



AvgMSE^slRp 



TMSE^B(Rp) 

iV2 



1 



1 



— + — )^ (24) 



(I24I 1 verifies the common sense that the more samples col- 
lected, the more accurate estimation acquired. And it also 
reveals that increasing SNR can reduce the estimation error. 
Furthermore, since 



— Xp rjXp — ||Xp||2 X 



X^Xp 



|xp||2 xR,Jf2) (25) 



where Rxp(r2) is the Rayleigh quotient of associated with 
the pilot sequence Xp, and Rx (11) < Xmax where X^ax is 
the maximum eigenvalue of ft. Besides, when the power of 
pilot symbol is normalized, l|xp||2 ~ N. Hence (l24l i is further 
lower bounded by 



AvgMSEi5(Rp) 



Nt 



(1 



(26) 



To further look into the relationship between f^Tg and Xmax, 
we examine the extreme eigenvalues of Jl for different fdTg's 
and A^'s numerically, and the results are plotted in FiglT] 
Moreover, we find a simple function fitting the maximum 
eigenvalues of all cases very well. The function is 



Xmaxi(n)=NJoi2TrcfdTs) 



(27) 
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Fig. 2. Comparison of analytic results 1241 and numerical results for EVA Fig. 3. Comparison of analytic results (28) and numerical results for EVA 
and ETU channels when 7 = 20dB and fa = 200Hz. and ETU channels when 7 = 20dB and fd = 200Hz. 



where c = 0.35 when fj^T^ < 0. Therefore, a more 
insightful lower bound can be achieved by 

1 1 



AvgMSE^s(Rp) = —(1 



N^M2TrcfdTsh 



(28) 



According to (1281 1, we know that the amount of samples, i.e., 
Nt, effects the estimation accuracy dominantly but SNR and 
maximum Doppler spread do not, since N"^ is sufficiently large 
for most of current systems. 

IV. Numerical Results 

The OFDM system in simulations is of BW = 1.25 MHz 
(T = 1/BW = 0.8 ms), N = 128, and L^p = 16. Two 3GPP 
E-UTRA channel models are adopted: Extended Vehicular A 
model (EVA) and Extended Typical Urban model (ETU) [18]. 
The excess tap delay of EVA is [0, 30, 150, 310, 370, 710, 
1090, 1730, 2510] ns, and its relative power is [0.0, -1.5, 
-1.4, -3.6, -0.6, -9.1, -7.0, -12.0, -16.9] dB. For ETU, 
they are [0, 50, 120, 200, 230, 500, 1600, 2300, 5000] ns 
and [-1.0, -1.0, -1.0, 0.0, 0.0, 0.0, -3.0, -5.0, -7.0] 
dB, respectively. The classic Doppler spectrum, i.e., Jakes' 
spectrum [11], is applied to generate the Rayleigh fading 
channel. 

In Fig|2] we compare the analytic results ( l24l i and the 
numerical results over a range of Nt's for EVA and ETU chan- 
nels, respectively, when 7 = 20dB and fd = 200Hz. The pilot 
sequences are QPSK modulated and randomly chosen. And 
the collected samples are apart from each others far enough 

*This condition ensures that Jo(27ra/^Ta) is positive and monotonically 
decreasing with respect to fdTs- In fact, this condition is always satisfied since 
current applied OFDM systems have f^Ts < 0.1 to maintain the power of 
ICI within a tolerable range [13]. 



to guarantee the assumption of independence. Apparently, the 
analytic results meet the numerical ones quite well. 

In Fig [3] we compare the analytic results ( |28T l and the nu- 
merical results for EVA and ETU channels, respectively, when 
7 = 20dB and fa = 200Hz. The pilot sequences are QPSK 
modulated. In order to examine the effect of different pilot 
sequences on w, one hundred different sequences randomly 
generated are tested and their MSE's are averaged and plotted. 
From the figure, we find that ( |28] l is a tight bound even for 
an arbitrary pilot sequence. 

The distributions of avgMSE for different SNR's and 
Doppler's are plotted in Fig|4] through ten thousands estima- 
tions for EVA and ETU channels, respectively. The amount 
of samples of each test is 200, and the pilot sequences are 
QPSK modulated and randomly generated. Clearly, avgMSE's 
are centered around zero and most of them are within the 
range of zero to CRLB, which follows that (fT6] l is an unbiased 
estimator Moreover, it is also obvious that the distributions of 
avgMSE for EVA and ETU channels are negligibly influenced 
by 7 and fd, which follows the analytic lower bound (|28l ). 

V. Conclusion 

In this paper, the maximum likelihood estimator and CRLB 
of the frequency correlation matrix for OFDM systems in 
doubly selective fading channels are derived and analyzed. 
Through the analyses, we obtain an insightful lower bound of 
average MSE, i.e., ( |28T l, and according to which, the amount 
of samples shows a dominant impact on the accuracy of 
estimation while SNR and maximum Doppler spread have 
relatively small effect when the number of subcarriers are 
sufficiently large, although increasing SNR and decreasing 
maximum Doppler spread can help to reduce MSE sUghtly. 
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Fig. 4. Distributions of average MSE for EVA and ETU channels when Nt = 200. 



References 

[1] O. Edfors, M. Sandell, J. van de Beek, S. Wilson, and P. Borjesson, 
"OFDM Channel Estimation by Singular Value Decompostion," IEEE 
Trans. Commun., vol. 46, pp. 931-939, July 1998. 

[2] Y. Li, L. Cimini, and N. Sollenberger, "Robust Channel Estimation for 
OFDM Systems with Rapid Dispersive Fading Channels," IEEE Trans. 
Commun., vol. 46, pp. 902-915, July 1998. 

[3] R Hoeher, S. Kaiser, and R Robertson, "Two-Dimensional Pilot-Symbol- 
Aided Channel Estimation by Wiener Filtering," in IEEE ICASSP 1997, 
April 1997, pp. 1845-1848. 

[4] B. Yang et al, "Channel Estimation for OFDM Transmission in Mul- 
tipath Fading Channels Based on Parametric Channel Modeling," IEEE 
Trans. Commun., vol. 49, pp. 467^79, March 2001. 

[5] M. Raghavendra et al, "Parametric Channel Estimation for Pseudo- 
Random Tile-Allocation in Uplink OFDMA," IEEE Trans. Signal Pro- 
cess., vol. 55, pp. 5370-5381, November 2007. 

[6] X. Zhao and T. Peng and W. Wang, "Parametric Channel Estimation by 
Exploiting Hopping Pilots in Uplink OFDMA," in IEEE PIMRC 2008, 
Cannes, France, September 2008. 

[7] D. Shiu, G. Foschini, and M. Gans, "Fading Correlation and Its Effect on 
the Capacity of Multielement Antenna Systems," IEEE Trans. Commun. , 
vol. 48, pp. 502-513, March 2000. 

[8] T. Wang et al, "Performance Degradation of OFDM Systems Due to 
Doppler Spreading," IEEE Trans. Wireless Commun., vol. 5, pp. 1422- 
1432, June 2006. 



[9] Y. Li and L. Cimini, "Bounds on the Interchannel Interference of OFDM 
in Time-Varying Impairments," IEEE Trans. Commun., vol. 49, pp. 401- 
404, March 2001. 

[10] X. Zhao and M. Yang and T. Peng and W. Wang, "An Analysis of 
the Bias-Propeily of the Sample Auto-CoiTelation Matrices of Doubly 
Selective Fading Channels for OFDM systems," Submitted to IEEE ICC 
09, Dresden, Germany, June 2009. 

[11] R. Steele, Mobile Radio Communications. IEEE Press, 1992. 

[12] R. Clarke, "A Statistical Theory of Mobile Radio Reception," Bell Syst. 
Tech. J., pp. 957-1000, July-Auguest 1968. 

[13] Y. Choi, P. Voltz, and F. Cassara, "On Channel Estimation and De- 
tection for Multicarrier Signals in Fast and Selective Rayleigh Fading 
Channels," IEEE Trans. Commun., vol. 49, pp. 1375-1387, August 2001. 

[14] T. Ratnarajah, R. Vaillancourt, and M. Alvo, "Complex Random Ma- 
trices and Rayleigh Channel Capacity," Commun. Inf. Svst., vol. 3, pp. 
119-138, October 2003. 

[15] D. Maiwald and D. Kimis, "Calculation of Moments of Complex Wishart 
and Complex Inverse Wishart Distiibuted Matrices," lEE Proc. -Radar 
Sonar Navig., vol. 147, pp. 162-168, August 2000. 

[16] K. Mardia, J. Kent, and J. Bibby, Multivariate Analysis. Academic 
Press, 1979. 

[17] G. Golub and C. V. Loan, Matrix Computations, 3rd ed. New York: 

Johns Hopkins University Press, 1996. 
[18] "3GPP TS 36.101 v8.2.0 - Evolved Universal Terrestrial Radio Access 

(E-UTRA); User Equipment (UE) Radio Transmission and Reception 

(Relasase 8)," 3GPP May 2008. 



